Genetic enhancement of Trichoderma asperellum biocontrol potentials and carbendazim tolerance for chickpea dry root rot disease management

Advances in biocontrol potentials and fungicide resistance are highly desirable for Trichoderma. Thus, it is profitable to use mutagenic agents to develop superior strains with enhanced biocontrol properties and fungicide tolerance in Trichoderma. This study investigates the N-methyl-n-nitro-N-nitrosoguanidine (NTG) (100 mg/L) induced mutants of Trichoderma asperellum. Six NTG (3 each from 1st & 2nd round) induced mutants were developed and evaluated their biocontrol activities and carbendazim tolerance. Among the mutant N2-3, N2-1, N1 and N2-2 gave the best antagonistic and volatile metabolite activities on inhibition of chickpea F. oxysporum f. sp. ciceri, B. cinerea and R. bataticola mycelium under in vitro condition. Mutant N2-2 (5626.40 μg/ml) showed the highest EC50 value against carbendazim followed by N2-3 (206.36 μg/ml) and N2-1 (16.41 μg/ml); and succeeded to sporulate even at 2000 μg/ml of carbendazim. The biocontrol activity of N2-2 and N2 with half-dose of carbendazim was evaluated on chickpea dry root rot under controlled environment. Disease reduction and progress of the dry root rot was extremely low in T7 (N2-2 + with half-dose of carbendazim) treatment. Further, carbendazim resistant mutants demonstrated mutation in tub2 gene of β-tubulin family which was suggested through the 37 and 183 residue changes in the superimposed protein structures encoded by tub2 gene in N2 and N2-2 with WT respectively. This study conclusively implies that the enhanced carbendazim tolerance in N2-2 mutant did not affect the mycoparasitism and plant growth activity of Trichoderma. These mutants were as good as the wild-type with respect to all inherent attributes.

Introduction Trichoderma spp. are among the most promising biocontrol agents used against numerous plant pathogenic fungi [1]. Trichoderma is known to stimulate plant health by protecting the plant via mycoparasitism, anti-microbial secondary metabolite production, as well as local and induced systemic resistance against invading pathogens [2]. There are several reports on Trichoderma species as biocontrol agent, however, only a handful of species like T. harzianum, T. afroharzianum, T. viride, T. asperellum, T. koningiopsis, and T. virens based formulations are used to manage the plant diseases [1,3].
Often Trichoderma species exhibits greater variability in antagonistic capacity and bio-stimulant action on plant pathogens and plants [4]. Some strains are more suitable for biological control of disease and others for stimulating crop growth and nutrient uptake [5]. The limited efficacy and inconsistent performance against targeted pathogens remains a major drawback in the field conditions [6]. Various strategies have been employed to enhance the efficacy and consistency of biological control activity. Recurrent selection of Trichoderma strains, and random mutations with UV light or chemical mutagenesis followed by subsequent selection has been a successful strategy to enhance efficacy and fungicide tolerance of Trichoderma species [7].
There have been several successful reports regarding the development of novel Trichoderma strains by using mutagenesis which has led to enhanced fungicide tolerance and biocontrol potential [1]. Pesticide tolerance in Trichoderma would be a prerequisite to reduce the fungicide resistance in pathogens and to promote the plant health [8]. Pesticide-resistant Trichoderma in combination with reduced fungicide application would help in reducing the detrimental effects of the latter on soil along with ensuring satisfactory levels of crop protection [9]. Moreover, reduced amount of chemical usage could weaken pathogen propagules for subsequent attack and give additive/synergistic effect on plant growth [10].
Among fungicides, carbendazim is one of the most extensively used benzimidazole (methyl benzimadazol-2-yl carbamate) for control of soil-borne diseases through seed treatment and foliar application [8]. Carbendazim resistance in pathogen populations have been detected in the field shortly after an intensive and exclusive exposure [11]. However, reversion of resistant populations has not been reported even after the selection pressure was removed [12]. Primarily, benzimidazole functions by inhibiting the polymerization of tubulin monomers into functional microtubules in mitosis and through sequence rearrangement of 5'-flanking region of β-tubulin gene in some species [13,14]. The benzimidazole tolerance in Trichoderma strains are yet to be ascertained for enhanced fungicide tolerance and large-scale field applications. There have been reports to confirm benzimidazole tolerance in Trichoderma with substitution of amino acids as a result of single and double mutation in β-tubulin tub2 region [8,15].
The aim of this study was to develop T. asperellum with enhanced biocontrol activity and carbendazim tolerant mutants of the species by induced mutation followed by characterization of biocontrol potential of mutants. The study also demonstrates dry root rot (DRR) disease control ability with half-dose of carbendazim and its potential in integrated management of soil-borne diseases under field condition.
The cultures were reactivated in potato dextrose agar (PDA) (HiMedia M096) medium and then grown in test-tubes containing 10 ml of the same medium. These cultures were stored in a refrigerator at 4˚C for future use in assays. For each experiment, samples taken from these test-tubes were cultured for five days in Petri dishes (90×15 mm) containing 20 ml of PDA medium.

Chemical induction and isolation of carbendazim tolerant mutants
Conidial suspensions of T. asperellum were obtained from sporulating 3 days old culture by rinsing with sterile distilled water. After centrifugation at 5000 rpm for 10min, concentrated suspension of conidia was added to N-methyl-n-nitro-N-nitrosoguanidine (NTG) (100 mg/L in 0.05 M Tris buffer pH 6) for 1 h in a shaker (250 rpm/10m) at room temperature. The NTG treated conidial suspensions were washed thrice in sterile distilled water by centrifugation at 5000 rpm for 10min. Conidial suspension (0.1 ml) were plated on PDA containing 10 μg/ml carbendazim (97%, Sigma Aldrich), 500 mg/L streptomycin and 150 mg/L Rose Bengal. The plates were incubated at 25±1˚C with 12 h light/dark photoperiod and growth of the colonies were observed daily. After 5 days, colonies which were tolerating 10 μg/ml carbendazim (97%) were isolated and sub-cultured for 10 times on PDA without carbendazim to test their stability [19]. Stable carbendazim tolerant mutants were selected and designated as NTG (N1, N2 and N3) mutants. Further, these mutants were tested at different concentrations of carbendazim (10, 20. . .120 μg/ml). Mutants which showed the maximum tolerance to carbendazim (i.e., N2) was selected as the parental strain for 2 nd round of NTG treatment by following the above protocol. After 2 nd round of mutation three more carbendazim tolerating mutants were picked up and mutants referred as N2-1, N2-2 and N2-3.

Screening and determination of sensitivity of the mutants
Dual culture assay. An actively growing mycelial disc (5 mm) of mutants and test pathogens (F. oxysporum f. sp. ciceri, R. bataticola and B. cinerea) were placed near the periphery of PDA amended Petri dishes and incubated at 25±1˚C with 12 h light/dark photoperiod. Inhibition of average radial growth of pathogen was calculated in relation to the growth on control [20]. Three replications were maintained for each treatment and Petri dishes containing test pathogen without Trichoderma served as control. Per cent inhibition of mycelial growth was measured by using below mentioned formula.
where I = Per cent inhibition of pathogen mycelium; C = Growth of pathogen in control; T = Growth of pathogen in Trichoderma inoculated plate.
Sandwich systems for volatile interactions. Sandwich plate technique as described by Li et al. [21] was used to test the volatile antifungal activity of mutants against the three chickpea pathogens. Actively growing 3 days old cultures (5 mm mycelial disc) of mutants and test pathogens were centrally inoculated in a separate PDA Petri dishes. Top of the Petri dishes were removed and base containing mutants and test pathogens were placed facing each other forming an atmosphere sharing set. Three replications were maintained for each treatment and incubated at 25±1˚C with 12 h light/dark photoperiod. Linear growth of test pathogens in inverted plates were measured when the test pathogens covered the Petri dish, and per cent mycelium inhibition was calculated by using above mentioned formula.
Maximum carbendazim tolerating mutants N2 and N2-2, along with wild-type (WT) were subjected for dry mycelia production at different concentration carbendazim. A 5 mm mycelial disc of these cultures were inoculated to potato dextrose broth (PDB) (HiMedia) containing carbendazim at 0, 250, 500, 750, 1000, 1250, 1500, 1750 and 2000 μg/ml. Flasks were incubated at 25±1˚C with 12 h light/dark photoperiod for 7 days without shaking. After 24 and 48 h of incubation, conidial germination and formation of mycelium in carbendazim amended flasks were observed under microscope (Olympus Corporation, Japan). After 7 days of incubation, cultures were filtered through Whatman (No.1) filter paper and mycelia was air-dried at 70˚C for 12h in hot-air oven. The dry weight (g) of the mycelium was recorded.

Greenhouse experiments for management of chickpea dry root rot disease
Development of talc-based formulations of mutants. Based on carbendazim tolerance, mutants N2, N2-2 along with WT was selected for chickpea dry root rot disease management against moderately resistant Annigeri and JG 62 cultivars under greenhouse conditions. Actively growing 5 mm mycelial disc of these cultures were inoculated separately into 250 ml conical flask containing 100 ml PDB and incubated at 25±1˚C with 12 h light/dark photoperiod for 7 days without shaking. Trichoderma conidia load was fixed to 10 8 /ml and mixed with autoclaved talcum powder at 1:3 proportion (S1 Fig). Formulations were air-dried at room temperature and colony forming unit (CFU) was calculated by following serial dilution technique on TSM (Trichoderma-selective agar medium) containing 100 mg/L Rose Bengal.
Bio-efficacy of mutants against chickpea dry root rot disease. Disease suppression ability of mutants N2, N2-2 and WT was evaluated against chickpea dry root rot disease in completely randomized block design (CRBD). Dry root rot susceptible pigeonpea cultivar Annigeri and JG 62 seeds were treated with Trichoderma formulations at the rate of 5 g/kg seeds and half-recommended dose (RD) of carbendazim (Bavistin-WP 50) i.e. at the rate of 1 g/kg seeds. Treated seeds were shade dried before sowing. Chickpea dry root rot sick pots were prepared by following the protocol mentioned in ICRISAT Chickpea Diseases Phenotyping Bulletin (S2 Fig). A total of eight treatments viz. T 1 (sick pot); T 2 (WT); T 3 (N2); T 4 (N2-2); T 5 (RD of carbendazim); T 6 (N2+0.5RD carbendazim); T 7 (N2-2+0.5RD carbendazim) and T 8 (Healthy pot) were levied separately on Annigeri and JG 62 cultivar. Seven seeds were sown per pot (13x8x9 cm) in 3 replications. Temperature (35˚C) of the greenhouse and moisture levels of the pots (60%) were carefully maintained as per the procedure given in the Bulletin. The percent disease incidence (PDI) was recorded at every four-day interval from the day of first disease appearance till complete wilting of all plants in the pots. The percentage of diseases incidence was calculated using the following formula.

Per cent disease incidence ¼ No:of wilt infected plants Total of plants assessed X 100
The dry root rot disease progress curves were constructed from the PDI data and the areas under the disease progress curve (AUDPC) and apparent infection rate (r) were calculated [23].
Apparent infection rate was estimated by regressing the logit of disease proportion (decimal from of per cent disease severity) on time in days [24].
Assessment of translational changes in β-tubulin gene of mutants DNA extraction and amplifications of β-tubulin gene. β-tubulin gene (tub1 and tub2) was amplified from mutants N2, N2-2 and WT (S1 Table). The total genomic DNA was extracted from actively growing 3-days old mycelium using CTAB Method (Mukherjee et al. 2019) tub1 and tub2 genes were amplified in a 25 μl PCR reaction for 35 cycles (S2 Table). Amplified PCR product was electrophoresed in 1.5% agarose gel and bands were visualized under UV transilluminator.
PCR clean-up and gel extraction were followed by using Macherey-Nagel gel elution kit and was ligated with pJET1.2 blunt vector (Thermo Scientific Clone pJET PCR Cloning Kit). The ligated product was used to transform competent E. Coli DH5-α. Positive colonies were confirmed by using colony PCR. The plasmid was isolated using Macherey-Nagel plasmid DNA purification Kit, and vacuum air-dried. The clones were checked by Sanger sequencing (Macrogen Inc., South Korea).

Sequencing and homology modelling of β-tubulin gene
All the obtained sequences were aligned and checked for sequence variation using Clustal omega (https://www.ebi.ac.uk/Tools/msa/clustalo/) and the alignment was visualised by using ESPRIPT 3.0 software (https://espript.ibcp.fr/ESPript/ESPript/) [25]. The β2-tubulin amino acid sequences were retrieved from the NCBI-Protein database in FASTA format and used as a reference for the alignment. The alignment was then analysed for occurrence of deleterious mutations, or amino acid substitution by using Provean (http://provean.jcvi.org/seq_submit.php) online tool. Homology modelling was carried out with the SWISS-Model interface with the β-tubulin X-ray crystallographic structure (PDB: 4ihj) as a template. All the models were evaluated by the QMEAN4 score. Global model quality estimation (Z score −1.22 and the QMEAN score 0.692) suggested that the modelled protein is of good quality. The homology models of mutated β-tubulin were built by altering the appropriate parent amino acids using the WHAT IF program (http://swift. cmbi.ru.nl/servers/html/index). The mutated models were subjected to overall energy minimization by steepest descent method, using GROMOS 43B1 force field. PyMOL (http://www.pymol. org/) was used for inspection of models and crystal structures while CastP server (http://sts. bioengr.uic.edu/castp/calculation.php) was used for pocket calculations using probe radius 1.4 Å.

Statistical analysis
The data of antagonistic, volatile inhibition and abiotic stress (NaCl salt and temperature) tolerance of mutants were arcsine transformed to make residuals normal and then back-transformed for the presentation of the results [26]. The transformed data was subjected for analysis of variance (ANOVA) to know the level of significance at p<0.01 by using PROC GLM, SAS 9.4 (Statistical Analysis Systems Institute Inc. 2016). Significance of mean differences within mutants was tested by Duncan multiple range test at p<0.01 level of probability. Boxplots were generated to visualize the per cent inhibition of chickpea pathogens, salt and temperature tolerance in mutants using R statistical program (R Development Core Team 2020). The EC 50 value of carbendazim tolerance of mutants were analysed by using "DRC" and "MULTCOMP" library in R statistical program.

Evaluation of altered morphology of mutants
The NTG induced mutants were significantly differed from the wild-type in growth rate and most of the mutants covered the complete Petri plate within 96 h of incubation. Highest average mycelial growth was observed in WT (55.70mm) followed by N2-3 (54.40mm) and lowest growth was observed in N1 (49.23mm) and N3 (49.478mm) mutant respectively (Table 1). Mutants did not show any variation in morphology, sporulation, colony pigmentation, conidia and chlamydospores production level when compared to wild type.

Evaluation of carbendazim tolerance and mycelial biomass production in mutants
After 2 nd round of selection, it was noticed that mutants had surpassed natural tolerance to field dose of carbendazim (Bavistin-WP 50). The minimum inhibitory concentration (MIC) of the fungicide was calculated to determine the mutants' tolerance capacity. When the carbendazim tolerant mutants were grown at graded levels of carbendazim, mutant N2-2 and N2-3 was observed to be more compatible at high concentration of carbendazim (Fig 1). The 6 NTG (3 each from 1 st and 2 nd round selection) treated mutants were developed. The mutants which were resistant to 10 μg/ml of carbendazim were tested at graded (10-120 μg/ ml) levels of carbendazim. Mutant N2 showed maximum carbendazim tolerance compared to N3 and N1 mutants (Fig 2A). Further, N2 mutant was subjected for 2 nd round of NTG irradiation and 3 more (N2-1, N2-2 and N2-3) mutants were recovered. At higher concentration, the growth of mycelium was negligible or almost nil in N2-1 as compared to N2-2 and N2-3 mutants (Fig 2B). The effective concentration (EC 50 ) values of these mutants were calculated where among the mutants, N2-2 (5626.40 μg/ml) showed the highest EC 50 values followed by N2-3 (206.36 μg/ml) and N2-1 (16.41 μg/ml).
Mycelial dry biomass production was significant at p<0.01 per cent in the mutants against different concentration of carbendazim. Maximum dry biomass production was noticed in N2-2 and results were at par with carbendazim at 0, 250, 500, 750, 1000 and 1250 μg/ml ( Fig 2C). As the fungicide concentration increased, the mycelium biomass and conidial germination was also reduced at different concentration of carbendazim (S4 Fig).

Disease control potentials of mutants
The disease control potential of previously observed carbendazim tolerating N2 and N2-2 mutants from 1 st and 2 nd round of selections were evaluated under greenhouse conditions with controls. The results obtained in the two separate experiments are shown in S3 Table. Mutant N2-2 along with half-dose carbendazim treatment (T 7 ) in both the cultivars, were registered with less disease incidence and higher plant growth compared to other treatments against chickpea dry root rot. Even after 60 days of sowing, adherence of Trichoderma mycelium was observed in the chickpea roots of T 7 treatment (Fig 3). Treatment T 7 (N2-2+0.5 RD carbendazim) showed significantly lower AUDPC incidence (1980.95 and 1895.24) and apparent infection rate (0.11 and 0.11) in both the cultivars. Treatment T 3 (N2) and T 4 (N2-2) without half-dose carbendazim, showed the maximum AUDPC of incidence and apparent infection rate. There were no significant differences in case of disease progress and incidence between the cultivars Annigeri and JG 62 (Table 3).

Sequence variations in the β-tubulin genes of carbendazim tolerance mutants
N2 and N2-2 mutants were analysed to determine the sites of mutation in β-tubulin genes. Alignment of a nucleic acid sequences of tub1 and tub2 genes in WT, N2 and N2-2 revealed that the mutations occurred only in tub2 (protein length: 563 amino acid), whereas no mutations was found in tub1 gene. The amino acid sequence demonstrated that there were over 211 observed mutations (substitutions) in N2 and 275 in N2-2 out of which 17 and 14 were found to be deleterious, respectively (S4 Table). Substitutions occurring at positions 156,157,170,171,189,190,237,238,239,247,252,253,255,277,278,279 and 362 in N2 were found to cause deleterious mutations which implied that the mutations were non-synonymous in nature.   (Fig 4). This might be the result of synonymous/nonsynonymous substitution of amino acids which would be responsible for mutation and ultimately carbendazim resistance.

Discussion
Biological control with improved Trichoderma strains has proven to be a potential alternative to chemicals of several soil-borne plant pathogens [27]. In many cases, individual use of either  biocontrol agents or chemicals could not provide satisfactory levels of control over soil-borne [8]. Attempts have been made to obtain resilient fungicide resistant strains for usage in integrated disease management [8]. The chemical irradiations followed by assessment of mutant's biocontrol efficacy have been guiding the selection of superior mutants for field application [1,28]. In the light of enhancing biocontrol efficacy in T. asperellum, NTG induced mutants confirmed the fungicide tolerance which were morphologically similar to the wild-type [1]. In chemical mutagenesis, the recommended dose of mutagenic compounds should be as large as possible with minimum deleterious effects [29]. In our study, the optimal dose of NTG (100 mg/L) was sufficient enough to cause the mutations leading to carbendazim resistance. Repeated NTG induced mutation in carbendazim resistant mutant further increased the carbendazim tolerance. Possible reason for the difference in fungicide resistant among repeated mutants could be due to acclimatization to the shock events. Production of cell wall degrading enzymes (CWDEs) and consortium of antifungal compounds are known to be the most important components of biocontrol by Trichoderma spp. [30]. Induced mutations may or may not increase the mycoparasitism activity of the Trichoderma spp. [31]. We found that mutants were capable of inhibiting mycelium of F. oxysporum f. sp. ciceri followed by B. cinerea and R. bataticola as compared to wild-type. This is in agreement with Abbasi et al. [28] who reported that Trichoderma mutants can increase inhibition of hyphal growth of F. oxysporum, R. solani, R. bataticola, P. ultimum, M. phaseolina, S. sclerotiorum and F. grminearum mycelium in both dual culture assays.
It was demonstrated that Trichoderma spp. usually use multiple biocontrol mechanisms, rather than a single mechanism to suppress plant pathogens [30]. The arsenal of volatile metabolites produced in Trichoderma spp. have shown its ability to inhibit the plant pathogens and promote plant growth [32]. We found that volatile metabolites produced by mutants were effectively inhibited the mycelial growth of F. oxysporum f. sp. ciceri followed by B. cinerea and R. bataticola. Reduction of mycelium growth in these pathogens are in accordance with Faheem et al. [33]. The study reported that volatile metabolites from T. viride (Tv-1) were effective against F. oxysporum, R. solani, S. rolfsii, S. sclerotiorum, C. capsici, H. oryzae, A. brassicicola mycelium inhibition. It has been observed that Trichoderma spp. may produce specific volatile metabolites for each type of antagonized pathogen.
Quite often biocontrol agents fail to perform as expected against targeted crop disease. Combination of different disease management methods have significant impact on many perspectives. From practical point of view, combined application of biocontrol agents and halfdose of recommended fungicides give additive or synergistic effects on plant diseases [10] by weakening the pathogen propagules [34]. Carbendazim resistance in NTG induced mutants are stable. We found that linear growth and mycelial biomass production was more and highly compatible in N2-2 mutant. T. harzianum mutant developed by Viji et al. [35] showed carbendazim tolerance up to 30 μg/ml. Similarly, NTG induced T. viride mutants exhibited benomyl tolerance up to 50 μg/ml [19] and Jayaraj and Radhakrishnan [36] demonstrated 100 μg/ml tolerance in T. harzianum. Further, Mech [37] reported that mutants of T. harzianum can grow even at 0.1% carbendazim.
Combination of N2-2 with half-dose of carbendazim treatment (T 7 ) significantly reduced the dry root rot disease incidence in Annigeri and JG 62 cultivar. Since, these genotypes are known to be moderately resistant, the T 7 treatment has effectively enhanced the resistance capacity against DRR. JG62 is known to exhibit a mortality rate of 90% under DRR incidence placing them sometimes under highly susceptible genotypes as well [38]. AUDPC and apparent infection rate (r) indicated slow disease progress in this treatment. R. bataticola in legumes is becoming intense in tropical humid areas [39] moisture deficit conditions and higher soil temperatures which are predisposing factor for disease in chickpea [40].Moreover, carbendazim tolerating T. harzianum mutants have been found superior in enhancing the root biomass of chickpea and increasing F. oxysporum f. sp. ciceri management to about four to eight fold [41].
In this study, we demonstrated that only tub2 gene shows mutations in N2-2 and N2. The deduced amino acids structural superimposition of tub2 gene in N2 with WT and N2-2 with WT, differ by 37 and 183 mutant residues. Carbendazim is known to inhibit the synthesis of β-tubulin protein [42]. Mutations in the β-tubulin gene evade the inhibitory effect of carbendazim and confer carbendazim tolerance to the mutant. A single amino acid change from phenylalanine to tyrosine at position 167 of the β-tubulin gene results in benomyl tolerance in N. crassa and A. chrysogenum [43]. Goldman et al. [15] confirmed the MBC resistant tub2 gene from T. viride strain T9BR47. A single mutation at amino acid 168, having Phe (TTC) instead of Ser (TCC) and a double mutation in amino acid 13 resulting in the substitution of Gly (GGC) by Val (GTG) in tub2 confirms the resistance in T. harzanum. Previous studies have demonstrated the mutations occurring at multiple codons of the β-tubulin gene ultimately confer tolerance to benzimidazoles such as carbendazim [8]. Our study reported 17 mutations in N2 and 14 mutations in N2-2 which were non-synonymous in nature.
This suggests that carbendazim tolerance could be attributed to mutations in tub2 gene in N2 and N2-2. Although N2 and N2-2 mutants showed increased carbendazim resistance, their biocontrol properties that prevented dry root rot were not hampered. This suggests that the fungicidal and biocontrol characteristics of Trichoderma are not linked to naturally occurring tub2 gene. This further shows that NTG induced mutagenesis allowed for selection of carbendazim tolerant Trichoderma mutants for management of dry root rot in chickpea. Further structural analysis of mutant tub2 protein structures would be necessary to narrow down upon the causal mutation.
Our results showed that the fungicide tolerant Trichoderma strain with half-dose fungicide could control major chickpea fungal pathogens without affecting the original essential mycoparasitism and plant growth activity. These mutants were as good as the wild-type with respect to all of these inherent attributes.

Conclusions
In conclusion, optimised mutagenic techniques improves the biocontrol potentials of Trichoderma. This study provides insight into the biocontrol potential of T. asperellum and its carbendazim tolerance. Additionally, the study demonstrates the soil-borne disease control ability of Trichoderma with half-dose of fungicide and confirms the carbendazim resistant mechanism in N2-2 mutants. The fungicide resistant Trichoderma mutants would certainly compliment the disease management practices with great environmental relevance.
To elucidate the antagonism and disease control potential of mutants with half-dose of fungicides against major chickpea pathogens, its molecular and biochemical basis of resistance needs to be examined. Whole genome sequences could explain the possible genetic causes for the heritably-stable changes in our Trichoderma mutant's antagonism and fungicide tolerance. There could be single and multiple point mutation for enhancing cell wall degrading enzymes and volatile metabolites productions and carbendazim-fungicide resistance in the mutants. Furthermore, experiments on tripartite interaction system that includes plant, plant pathogen and biocontrol agent (consortium of improved antagonistic and fungicide tolerant mutants) would provide the significance of improved biocontrol agents specifically in additive or synergistic effects on plant diseases, plant growth promotion, soil health as well as lowering chemical residue in crops and soil.